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1. Introduction 

In this paper, we consider a lattice consisting of M spin 5 = 1, interacting with 
each other through a pairwise interaction. An example is the Hamiltonian for a one- 
dimensional spin-1 chain with nearest neighbor interaction, 



where the sum is over the site labels I and the second relation defines 7. This 
Hamiltonian has been of intense interest in the theoretical studies on quantum 
magnetism PI El El E] Here we shall demonstrate some exact properties of the ground 
states for this type of Hamiltonians. For definiteness, we shall confine ourselves to the 
Hamiltonian ([TJ (or its slight generalization (JHJ) below) with M being even and finite 
unless otherwise stated. We shall be particularly interested in the region K < J < 0, 



In a recent paper we have pointed out that this region of 7 is relevant to 
spin-1 Bosons trapped in an (optical) lattice in the regime of one particle, one orbital 
state per site, for suitable interactions between the Bosons. To see this, we generalize 
the standard derivation of the Heisenberg (exchange) Hamiltonian from the Hubbard 
model P to spin 1. Let the hopping between sites be t and the (repulsive) interaction 
for two particles on the same site be U (> 0) if their total spin is 0, and U 2 (> 0) if 
their total spin is 2. (Note that they cannot have total spin 1 since we have identical 
Bosons and there is only one orbital per site.) Consider now two neighboring sites (say 
1 and 2) with one Boson each. At t = the states corresponding to different spin 
configurations are degenerate. For finite but small t we can perform perturbation to 
second order. The energies of the system can be classified according to the total spin 
Stat f° r the two sites, and can be easily seen to be — 4t 2 /[/o, 0, and —At 2 /U 2 for total 
spin S = 0,1,2 respectively. Comparing these values with those of the general form of 
the spin Hamiltonian e + J (Si ■ S 2 ) + K(Si ■ S 2 ) 2 for these two sites, and noting that 
Si-S 2 = -2, -1, 1 for St* = 0, 1, 2 respectively, we find J = -fg, K = -§£ - f £ 
and e = J — K. Thus f^<J<OifO<£/o< U%. (e$ represents an energy shift 
independent of the spin configurations and can thus be dropped from our discussions.) 
The case < Uq < U2 is applicable to 23 Na atoms. [7] 

Our major results are collected as Theorems and Corollaries in Sec 121 In particular, 
we shall show that the ground state for the Hamiltonian ([TJ has zero total spin S tot = 0. 
We also show that the excitation of lowest energy but with finite spin has >Stot — 2. 
Generalizations of our results to lattices other than one-dimensional chains are also 
possible. Some of them are mentioned in subsection 12.21 

For comparison, we shall, in section discuss the case where J < K < 
[-7T < 7 < This case applies to spin-1 Bosons in a lattice if the interaction 



n = Y, [J (Si- s z+ i) + k (s t ■ s m ) 2 



(i) 



= V J 2 + £ [ cos 7 (S l ■ S l+ i) + sin 7 (S, • S m ) 2 



(2) 



i.e., -f <7<-|- 
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among Bosons is "ferromagnetic", (L^ < Uq) as would be the case for 87 Rb atoms 
trapped in an optical lattice It is worth remarking here that spin-polarized 87 Rb in 
an optical lattice in the (Mott) regime of one particle per site has already been obtained 
in a recent experiment 

Our work generalizes the corresponding results for the Heisenberg ( J ^ 0, K = 
) Hamiltonian. For a bipartite lattice and J > (the antiferromagnetic Heisenberg 
model), the results are usually referred to as the Marshall theorems jHj- The proofs here 
follow similar ideas, but the details are very different (in particular, our results here do 
not rely on a bipartite lattice). Section 0] below contains a comparison between the two 
cases. Main results are summarized in section EJ 



2. The region K < J < 0: 

To ease the presentation we shall first consider a linear spin chain (subsection 12. lj) . Then 
we shall discuss more general lattice types (subsection 12. 2|) 



2.1. Linear spin chain 

We shall consider, for definiteness, a generalized form of eq (JTJ: 

H = E (l " [ J (Sj • S,+i) + K (S, ■ S m ) 2 ] (3) 

i 

where < 5 < 1, i.e., with bond alternation ( in strength but not in sign of the 
interaction). This is shown schematically in Fig^ For convenience we shall also define 
e = 1 — 5 > 0. eis thus proportional to the strength of the weaker bonds. For 5 = 
(e = 1), Hamiltonian (J3J) reduces to that of (0). For 5 = 1 (e = 0), the Hamiltonian can 
be solved trivially since we have a collection of spin pairs with no interaction among 
different pairs, a fact that we shall take advantage of. 

We shall use the 'Ising configurations' as our basis set. Each of these configurations 
is defined by specifying the spin projection along the z direction for every site. Denoting 
the state at a site by |+ >, |0 >, |— > according to S z = 1,0,-1, a typical Ising 

configuration, which will be represented by \a >, is, e.g., | H — OH 00... > etc. We 

shall express any many body state |^ > as a linear combination of these basis states, 
i.e 

l^>=E c «l a> ( 4 ) 

a 

where the sum is over all \a >. We shall also introduce the rotated basis, defined by, 
for all sites, 

1+ >' = 1+ > 

!->' = h> (5) 
10 >' = % 10 > 
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and also use the expansion 



(6) 



The factor % in eq © should not alarm us. Since the z component of the total spin 
M = Stot,z is conserved under TC, for each given M, the number of sites I with S Zt i = 
in the expansions (|1J) or (jOJ) must be either even or odd for all a's [ for given M, the 
number of these sites can be only changed by replacing H — pairs ( not necessarily 
nearest neighbours) by 00, or vice versa ]. 

Theorem 1: For the lowest energy state of TC of a given M, the coefficients c a > are non- 
negative for all \a >' (apart from an overall common phase factor, a caveat that we 
shall not repeat). 

Proof. The basic idea of the proof is similar to that of the Marshall theorem j^]. One 
regards the Schrodinger equation TC\^f >= E |\& > as the corresponding one for a 
tight-binding (hopping) Hamiltonian with the lattice points labelled by \ot > ( \a >' ). 
Consider the general term in the Hamiltonian between sites / and I + 1: 

J (Si ■ S m ) + K (Si ■ S m ) 2 



H u+l = (1 ± 5) x 



This term affects only the spins I and / + 1. When Hii + i operates on an Ising 
configuration \a >, spin configurations for V ^ 1,1 + 1 are unchanged. For simplicity, 
we shall suppress the configurations for all these /'. For the spin configurations of I and 
/ + 1, Hii + i only makes 'hopping' among \ -\ — >, |00>, | — h>; between | + > and 
1 04- >, and between | — > and |0— >. In the first mentioned subspace ( | 4 — >, 
1 00 >, | — h >), the matrix elements of -ff^+i are given by, 



(1±S) x 



V 



-J + 2K 
J-K 
K 



J-K K 

2K J-K 
J-K -J + 2K 



(7) 



In the second mentioned subspace 

' K J \ 
J K j 



4- > and |0+ > ), they are 



[l±S)x 



(8) 



with identical matrix elements for the third subspace ( | — > and |0- 

diagonal in | 4- 4- > and | >. 

In the rotated basis, the matrix ((7j) transforms to 

( -J + 2K -(J-K) K \ 



> 



1,1+1 is 



(1±5) x 



V 



(J-K) 
K 



(9) 



2K -(J-K) 
-(J-K) -J + 2K J 

while the matrix (JHJ) is unchanged. In our 7 region of interest, J — K > 0, K < 0, 
J < 0, hence all hopping (off-diagonal) matrix elements are < 0. It follows that the 
lowest energy state must have c' a > for all \a >'. 

Corollary 1: If 5 7^ 1, the lowest energy state for given M has c a i > for all a n s. It 
also follows that this state is unique for given M. 
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Proof: If 5 ^ 1 (e > 0), all hopping matrix elements are negative (non-zero). As in 
the case of the Heisenberg Hamiltonian [Bj, all Ising configurations are connected by 
(though multiple) hopping. By the same reasoning as in [Oj, the lowest energy state 
must have all coefficients d a non-zero; and moreover there can be only one such state. 
Remarks: For 5 = 1 (e = 0) the system becomes a collection of spin pairs, with no 
interaction among pairs. For one pair, the energies Eg only depends on the total spin 
S and are given by 



For K < J < 0, Eq < E2 < E\. The minimum energy state is given by, for 
M = 0, ^ [| + - > +| - + > -|00 >] = ^ [| + - >' +| - + >' +|00 >']; for M = 1, 
75 [I + > +|0+ >] = --*= [| + >' +|0+ >']. and for M = 2, | + + >= | + + >'. The 
first state has S = and the latter two both belong to S = 2. These states obviously 
obeys the signs stated in Theorem 1. 

For Af = 2N spins, the ground state has S t ot = M = and is a collection of singlet 
pairs. For M = 1 (2), the lowest energy states have one pair of spins in the S = 2, 
S z = 1 (2) state with the rest in the singlet states. These states are not unique due to 
the freedom of choice of which pair being the S = 2 pair. 

Theorem 2: For a given M, the lowest energy state for ^ 5 ^ 1 has the same Stot 
independent of 5. 

Proof: For any given 5 7^ 1 we have already seen that the minimum energy state has 
c a i(5) > for all Ising configurations a'. For 5 = 1, we have c a >(l) > 0. Their overlap, 
given by J2 a ' c «'(l) x c a'(8), is non-vanishing. Hence they must have the same S^t- 
From the Remarks following Corollary 1, it follows that, for our region of 7, 

(i) The lowest energy state in the M = sector has S tot = 0, 

(ii) the lowest energy state in the M — 1 sector and the lowest energy state in the 
M = 2 sector both have S^t = 2. 

Since the M sector contain states with S tot > M, we have 
Corollary 2: 

(a) The ground state for our region (K < J < 0) has S^t = 

(b) The spin excitation (5* 7^ 0) with the smallest excited energy belongs to Stot = 2. 

(c) Any state with S^ot = 1 has energy higher than the mentioned state in (b) if 5 ^ 1 



The lines of reasoning above can also applied to deduce some properties concerning 
states of large S'tot. The sector M = 2N — 1 has one state belonging to S^ot = 2N and 
one belonging to S'tot = 2N — 1. The former state can trivially be written down and 
has c a i > 0. Hence we conclude that E 2 n < E 2 n-i, and this later expression becomes 
a strict inequality if e 7^ 0. 

For odd Af, similar considerations above show that the ground state has S to t = 1- 
The conclusions in this Corollary agree with known numerical results (e.g. [3]). 



E 
E 2 



2 [-2 J + 4K] 
2 [- J + K\ 
2 [ J + K) 



(10) 



(e>0) 
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2.2. Other lattice types 

As may already have been obvious, our Theorems and Proofs in the last subsection can 
be generalized to other lattice types. Moreover, these lattices need not be bipartite. A 
particular interesting case is as shown in Fig^D. Here we are again considering a spin-1 
lattice with pair-wise interactions J (Si ■ S 2 ) + ^(Si ■ S2) 2 with the strengths of the 
bonds proportional to, for full lines, 1, dashed lines, ei, and dotted lines, e 2 . (ei,e 2 > 
but not necessarily equal.) The case where all bond strengths are equal is contained 
in ei = e 2 = 1. The arguments leading to Theorem 1 and Corollary 1 went through 
unchanged. If ei = e 2 = 0, the system again reduces to the collection of pairs, and the 
Remarks subsection following Corollary 1 applies. By comparing the system of interest 
(with €1, e 2 both finite) with that with e.\ = e 2 = 0, we can again prove Theorem 2 and 
Corollary 2 as before. It is worth mentioning that a triangular optical lattice can also 
be formed by three suitable laser beams propagating at angles 7r/3 with respect to each 
other, thus the discussions here is applicable also to a physical system. 

3. The region J < K < 

Similar arguments above can easily be generalized to the region J < K < (— 7r < 7 < 
— ^f) for Hamiltonian (J3J). In this case with < 8 < 1, all hopping matrix elements 
are negative in the original (unprimed) basis, and thus the ground state for each M 
sector has c a > 0. This result is in agreement with the common wisdom expressed in 
the literature (e.g., Ref 0j) that the ground state here has S tot = 2N. These states, 
obtainable from | + + + +...> by suitable number of lowering operators, indeed has 
c a > for all \a >'s. 

4. Comparison with Marshall theorems 

For the antiferromagnetic Heisenberg (J > 0, K = 0) model on a bipartite lattice, the 
Marshall theorems hold. The basis employed in these theorems is defined by, for the A 
sublattice, 



+ >M 


= |+ > 


> M 


= |0 > 


— >M 





and for the B sublattice, 

1+ >M = ~ 1+ > 

|0 > M = |0 > 

|- >M = - |- > 

For our more general Hamiltonian (j3J), 
is given also by (JOJ) for the | H — >, 1 00 >, | 



(12) 

the matrix elements of Hij + i in this basis 
— h > subspace while that for I + > and 
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|0+ > (or | — > and |0— >) becomes 

<- 1±s >*{-j~k) (13) 

Thus, for J > 0, K < (J — K > trivially) all "hopping" elements are negative. 
Thus the Marshall theorems can be generalized to the region — ~ < 7 < with the 
results (arguing as in |Hj or as in Sec EJ note that in this region, for a system with two 
spins, one has E < E\ < E 2 ) 

(&M) The ground state of (jSJ) has S tot = and is unique ( if 6 ^ 1) 
(bM.) The spin excitation with lowest energy has S to t = 1- 

(aM) Any state with S to t = 2 has energy higher than the state mentioned in (bAi) if 
e > 0. 

For 7 = — 7r, some hopping matrix elements vanish (J = 0). There are in general 
degeneracies for the minimum energy states for a given M. This fact can also be seen 
from the work of Parkinson j2] . 

5. Concluding remarks 

We have proven some exact properties of lowest energy states for Hamiltonian of the 
type (JHJ) with K < J < 0. We demonstrated that, for a large class of lattices (not 
necessarily bipartite) with even number of sites, the ground state belongs to total spin 
^tot = 0, whereas the state of minimum excited energy but with finite S'tot belongs to 
S'tot = 2. The results derived here are relevant to spin-1 Bosons trapped in an optical 
lattice in the regime of one particle per site for suitable interaction between the Bosons. 
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Figure 1. (a) Schematic representation of the spin 1 chain with Hamiltonian 
Thick (dotted) lines represent stronger (weaker) bonds, with strengths proportional to 
1 + 5(1—5). (b) Schematic representation of a hexagonal spin 1 lattice. 



